A Method for Modeling Decoherence on a Quantum Information Processor 
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We develop and implement a method for modeling decoherence processes on an N-dimensional 
quantum system that requires only an iV 2 -dmiensional quantum environment and random classical 
fields. This model offers the advantage that it may be implemented on small quantum information 
processors in order to explore the intermediate regime between semiclassical and fully quantum mod- 
els. We consider in particular o z a z system-environment couplings which induce coherence (phase) 
damping, though the model is directly extendable to other coupling Hamiltonians. Effective, irre- 
versible phase-damping of the system is obtained by applying an additional stochastic Hamiltonian 
on the environment alone, periodically redressing it and thereby irreversibliy randomizing the sys- 
tem phase information that has leaked into the environment as a result of the coupling. This model 
is exactly solvable in the case of phase-damping, and we use this solution to describe the model's 
behavior in some limiting cases. In the limit of small stochastic phase kicks the system's coherence 
decays exponentially at a rate which increases linearly with the kick frequency. In the case of strong 
kicks we observe an effective decoupling of the system from the environment. We present a detailed 
implementation of the method on an nuclear magnetic resonance quantum information processor. 
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I. INTRODUCTION 

As early as the 1930s von Neumann 0] recognized that 
quantum correlations are crucial to understanding the 
quantum measurement process. He considered measure- 
ment as a process that first required correlating the sys- 
tem with the quantum apparatus through a unitary, in- 
formation conserving, quantum evolution. To complete 
the measurement a mechanism was needed by which this 
pure, correlated state decayed into a mixture approxi- 
mately diagonal in the basis of observation. In recent 
decades, the process of decoherence, which explains the 
dynamical origin of the above decay, has been extensively 
studied 0, 0, M 13 ■ By employing an open-systems ap- 
proach, the effect of the interaction between the measure- 
ment apparatus and its environment was included ex- 
plicitly, and von Neumann's method was extended. The 
physical origin of the process of decoherence is very sim- 
ple: the quantum correlations between the apparatus and 
the environment that are established in the course of their 
interaction is responsible for the dynamical selection of 
a preffered set of states of the apparatus (the pointer 
states). The mechanisms of decoherence are now a sub- 
ject of great practical interest. Some of the recent work 
on decoherence includes the determination of emergent 
properties of pointer states 0, , efforts to design spe- 
cific pointer states by engineering the environment [8|, 
and identification of the time-scales of the decoherence 
process ||. 

One of the simplest and most illustrative models of de- 
coherence was originally suggested and studied by Zurek 
H . It consists of a two level system (a spin 1 /2 particle) 
coupled to n two level systems through a cr z a z type inter- 
action. With this model, in the large n limit, it is possi- 



ble to show that the correlations which arise between the 
system and the environment lead to the damping of the 
system coherence, encoded in the off-diagonal elements of 
the density matrix. In this work we will present results, 
both theoretical and experimental, for a two-level system 
that is coupled to a few other two-level systems, which 
shows that by manipulating the latter one can reproduce 
the essential features of Zurek's model. 

Interest in this and other decoherence models (for 
example a two-level system coupled to a boson bath 
|10LlllUl2j ). has grown over the last few years due to the 
development of quantum information processing (QIP). 
A major challenge in QIP is the preservation of quan- 
tum coherence in the face of constant perturbations by 
an environment. While one could try to isolate the QIP 
device, this would make controlling the system difficult. 
Therefore, other strategies like quantum error correction 
(QEC) [Hi and noiseless subsystems (NSs) [H E3 
have been developed. The aim of this work is to develop 
methods to emulate decoherence in a physical setting, 
such as a QIP device, so that the nature and underly- 
ing physics of decoherence can be better understood and 
applied in the development of control strategies. 

The paper is organized as follows. In Section II we 
introduce the essential features of decoherence reviewing 
the model proposed by Zurek in which the system con- 
sists of a single spin while the environment is composed 
of an ensemble of n spins. In Section III we describe 
a simple model in which the environment is limited to 
only a few spins (qubits) and analyze a strategy through 
which these few spins can simulate a much larger effec- 
tive environment. The strategy consists of randomly re- 
dressing the phase of the environment qubits during their 
interaction with the system and averaging over many re- 
alizations of this evolution. We describe an exact solu- 



2 



tion of this model in the case of a a z <j z coupling between 
the system and a single environment qubit. In this case 
we provide an analytic description of the decoherence 
(phase-damping and decoupling) effects that arise under 
specific limiting conditions and also derive the associated 
Kraus operators for the model. A more detailed numer- 
ical analysis of this model is given for the case in which 
the environment consists of two qubits. In section IV, 
we present nuclear magnetic resonance (NMR) QIP sim- 
ulations for the two qubit environment and comparisons 
of these results with the one and two qubit environment 
predictions and numerical simulations. In Section V, we 
summarize our results and discuss the extension of this 
model to more general decoherence mechanisms. 

II. ZUREK'S DECOHERENCE MODEL 

In this section we review the basic elements of quantum 
decoherence by presenting an open-system model due to 
Zurek []| which is simple enough to be solved analytically. 
In spite of its simplicity the model captures many of the 
elements of decoherence theories and sheds insight into 
the loss of coherence, the onset of irreversibility, and in 
particular, the role played by the size of the environment. 

Consider n two-level systems and focus on one system 
as the subsystem of interest. This subsystem interacts 
with the rest of the system through a bilinear interaction. 
The overall dynamics is described by 

n 

'Use = Jik<jl<Jz, (1) 

k=2 

where the system qubit is denoted by the superscript T'. 
This Hamiltonian is energy conserving and only causes 
phase damping. The prescription of the open-systems 
approach is to evolve the combined system and environ- 
ment, represented by the density matrix p SE (t), and then 
recover the system density matrix from a partial trace 
over the environment degrees of freedom: 

/ (()= T rE{ /- W} =(|»^W). (2) 

In Eq. [21 Poo (0 an( i Pii M represent the system popula- 
tion terms while pQi(t) = pfo(t) represents the system 
coherence term. If the coherence terms vanish, the pure 
state is turned into a mixture in the computational ba- 
sis (er 2 -basis), i.e. a "pointer basis" has been selected 
out by einselection. An important observation is that, in 
the absence of a self-Hamiltonian, the system's statonary 
states are selected out by the interaction Hamiltonian. In 
fact, since [a z ,TLtot] = 0, the interaction with the envi- 
ronment has two memory states |0)s, \l)s as eigenstates 
and the populations remain unchanged throughout the 
system's evolution. The coupling in Eq. \T\is therefore a 



purely phase damping mechanism and there is no energy 
exchange between system and environment. 

The combined system evolves by the unitary propaga- 
tor 

n 

UsE{t) = exp(-iHsEt) = exp(-i^2 J lk ala^t). (3) 

k=2 

Consider a factorizable initial state of the combined sys- 
tem: 

M0)) SE = \1>(0))s®\1>(0))e 

n 

= (o|0>i+6|l>i) IJ(a fc |0) fc + /3 fc |l)fc). (4) 

k=2 

The evolution is such that 

n 

Mt))sE = a\0) 1 l[e- u ^ t \cf ) ) k 

n 

+b\l)x n e***^)*, (5) 

k=2 

where \4>)k = &k\0)k + /3fe|l)fe- The interaction entangles 
the system states with the environment. In the language 
of QIP, the transformation Use generates a conditional 
phase gate between the system and its environment, con- 
ditioned on the system's state. After the interaction the 
state is 

n 

m))sE = o|Q)i Y[[a k e- u ^\0}k + /3 fe e a " 4 |l) fc ] 

fe=2 

n 

+b\l) 1 l[[ake u ^ t \0)k + (3ke- u ^\l) k } (6) 

k=2 

and reflects the fact that the system and environment 
states are not factorizable. The off-diagonal element of 
the system's reduced density matrix (system coherence) 
is 

P S oS) - t(0\Tr E {m))sEm)\sE}\lh, (7) 

so that 

p§ 1 (t) = ab*z(t) (8) 

where 

n 

m = Y[[\a k \ 2 e- 2U ^ + \P k \ 2 e 2U ^}. (9) 

k=2 

Recall that a and b are the coefficients of the initial pure 
state of the system. 

The time-dependence of z(t) contains the crucial infor- 
mation for understanding the behavior of the system co- 
herence. In particular, the magnitude of z(t) determines 
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the damping of the phase information originally con- 
tained in poi(O), with \z(t)\ — > reflecting non- unitary 
evolution and "irreversibility" . 

For a finite system, \z(t)\ is at worst quasi-periodic and 
one can always define a recurrence time, te- The exis- 
tence of such a recurrence time reflects the fact that the 
information loss is in principle recoverable. In the con- 
tinuum limit, n — > oo, z(t) is no longer quasi-periodic 
and te — > oo. The phase information is then unrecov- 
erably lost, displaced from the degrees of freedom of the 
system to the infintcly many degrees of freedom of the 
environment. 

To characterize the degree of decoherence one can con- 
sider the size of the fluctuations of z(t) around its time- 
averaged mean value < z(t) >= 0: 



< z 



(t)\ 2 > = lim ± f dt'W)? 

1 n 

« ^ill^ + d^l 2 -^! 2 ) 2 ]- 



fe=2 



(10) 



Thus, typical fluctuations vary as ^==== and the effec- 
tiveness of the decoherence mechanism in this model is 
determined by the dimension of the environment. 

To summarize, the key features of this model of deco- 
herence are: (1) the system of interest evolves through 
a direct entangling interaction with each two-level sys- 
tem in a very large environment, and at any time the 
(reduced) density matrix of the system is obtained from 
a trace over the environment degrees of freedom; (2) ex- 
pressed in the pointer basis of the system, which in this 
simplest case is the set of states that commute with the 
interaction Hamiltonian, the reduced density matrix be- 
comes approximately diagonal and the off-diagonal ele- 
ments exhibit coherence loss; (3) the fluctuations of the 
decoherence produced by this model, measured by the 
size of the system's off-diagonal elements, are controlled 
by the dimension of the environment's Hilbert space. 



III. A HIERARCHICAL DECOHERENCE MODEL 

We consider the problem of experimentally simulat- 
ing quantum decoherence in a physical setting in which 
limited quantum resources are available for modeling the 
quantum environment. By "simulating quantum deco- 
herence" we are referring not only to the challenge of 
implementing an arbitrary open-system trajectory on a 
QIP device, but also to the study of the decoherence 
processes that result from specific system-environment 
couplings (for example, derived from a model of some 
physical system of interest). 

As in Zurek's model the exclusive direct mechanism for 
system decoherence in our model is the coupling between 
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FIG. 1: The schematics on the left describe the models de- 
veloped in this paper. A single spin system, S, is coupled 
to one and two spin quantum environments, designated Ei. 
During the coupling the system phase information leaks into 
the environment. Since the spin environment is finite, in or- 
der to simulate the effects of a larger quantum environment 
(depicted on the right) a mechanism is needed by which the 
information stored in the available quantum environment can 
be effectively erased. We accomplish this by redressing the 
environment degrees-of-freedom with stochastic phase kicks. 
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FIG. 2: A schematic of the system coupled to a hierarchy 
of quantum environments, as resources permit, and the role 
of the classical stochastic field. Each environment is coupled 
only locally (to its neighbours in the hierarchy). It should be 
noted that the classical environment only interacts with the 
quantum environments, and does not interact directly with 
the quantum system. In this paper we consider only the case 
of one local quantum environment, as portrayed in Fig. 1. 



4 



the system and a local quantum environment through a 
fixed bilinear Hamiltonian. However our model of de- 
coherence has two distinct features from the model de- 
scribed above. The first difference is a constraint on the 
Hilbert space size derived from practical considerations: 
we allow the dimension of the Hilbert space for the local 
quantum environment to be no larger than iV 2 , where TV 
is the dimension of the Hilbert space of the system. In 
this way the quantum environment is the smallest size 
that will enable the implementation of an arbitary com- 
pletely positive map on the system through a unitary 
operator on the combined system and environment. To 
remove the information from our finite quantum environ- 
ment we include a stochastic classical field in our model. 
This strategy is designed to eliminate the quantum back- 
action from low dimensional environments. Basically the 
technique consists of redressing |l7| the environment's 
quantum state by applying a sequence of random clas- 
sical kicks to the environment qubits, and then averag- 
ing over realizations of this stochastic noise. This has 
the effect of scrambling the system information after it 
has been stored in the quantum environment through the 
coupling interaction. It is worth stressing that the sys- 
tem itself is not subjected to these classical kicks and 
the associated stochastic averaging. This model, and 
the associated method realized in this paper, is depicted 
schematically in Fig. 1. A generalization of this method 
to provide a time-dependent open-system evolution, is 
described in the discussion and depicted in Fig. 2. 

As we shall show below, this scheme enables simula- 
tion of the quantum decoherence that normally arises for 
much larger effective environment sizes. In particular, 
we demonstrate the simulation of phase-damping on an 
NMR QIP consisting of three qubits (see Fig. 1). In 
the NMR simulation, the system is represented by one 
qubit while the other two qubits represent the quantum 
environment. Before turning to a discussion of the three- 
qubit experiment, we first describe and analyze this sim- 
ulation method theoretically in the simplest and solvable 
case of the phase-damping of a single system qubit from 
a single environment qubit (also depicted in Fig. 1). 



One-qubit Environment: Simple Solvable Model 

Below we introduce the essential features of this deco- 
herence model by considering an exact solution available 
in the case of a one-qubit environment coupled to the 
system by a a z a z interaction. With the system and en- 
vironment qubits labeled by S and E, respectively, the 
full Hamiltonian is given by 



Wo = k(v S (Jz + v E a l z 



2 a z a z ) 



(11) 



Here, vs, v e and f2 are frequencies in units of Hertz. This 
Hamiltonian includes both the self-evolution of the two 



qubits and their interaction. In the absence of any other 
interaction, the evolution operator for a time t is 



U(t) = exp 



-m{vscr z + v E o\ 



9. 



■ S „E 



(12) 



We will consider the evolution of this system subject to a 
sequence of kicks that affect only the environment qubit. 
Every kick is generated by a transverse magnetic field 
that rotates the environment qubit around the y-axis by 
an angle e m chosen randomly in the interval (—a,+a). 
The evolution operator for the m-th kick is given by K\ = 
I s ® Kf, where 



Ki 



exp(- 



i£ m cr„ 



(13) 



and / is the identity matrix. In our proposed model, the 
kicks are considered instantaneous, therefore the evolu- 
tion for a total time T = n/T, where T is the kick rate, 
can be written, 



U n (T) = K n U [-) K n -!U (- 
n ] \n 



T 



(14) 



It is useful to keep in mind that the operator U n {T) 
depends also on the values of the random variables e m 
(to = 1, n) corresponding to the kick angles. 

Our goal is to obtain a closed expression for the re- 
duced density matrix of the system qubit for an ensemble 
of realizations of the random variables e m . The density 
matrix for this ensemble is given by 

^/It-'Xf^^ ^' (15) 

where a is the spread of allowed kick angles over which 
the e m (to = 1, ...,n) are uniformly distributed. We will 
consider a factorizable initial state for the two qubits 
(this is not essential): 



^*(0)=^(0)®p*(0). 



(16) 



It is convenient to express the initial density matrix of 
the system in the basis of eigenstates of a z , 



3,1=0,1 



(17) 



Then we can simplify the expression for U n (T). To do 
this, we evaluate the effect of the first step in the evolu- 
tion Eq. ^]as follows, 

p SE (l) = KxUiT/n) p SE (0) U{T/n)^K\ (18) 

= E [p»b'X'G 

® [K E V E p E (0)(K E Vf)t] , 
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where we have defined the environment operator, 



(19) 



exp 



-i-v s (-iy-i- U{-iy + v E )a E 



In the above we have explicitly evaluated the action of 
the interaction Hamiltonian on the system states, and 
the j-dependence of the single-step operator V E reflects 
the fact that it operates on the environment state con- 
ditionally on the system state. The important point is 
that the evolution operators for the additional n — 1 iter- 
ations will factor as above, producing a final expression 
with (conditional) operators that act exclusively on the 
environment qubit. Hence, we can immediately obtain 
the following simple form for the final density matrix of 
the system qubit, 



P S {T) 



Pjl(Q)fjl(n,T) 

3,1=0,1 



(20) 



where the function fji(n,T), which we call the decoher- 
ence factor, carries all the information about the effect 
of the environment qubit on the system qubit, including 
also the trivial phases from the system's self-evolution. 
It is given by the formula 



fji{n,T) 



f a de n 


r dei 




■J- a ^ TrE 


l- a 2a ' ' 



(Af) nP E (0) (Af)l' 
(21) 



= K E V E K E ,V E 



■ KtVf 



where the operator (A E ) n n is defined as 

(Af) n = (j\U n (T)\j 

" (22) 

It is clear from Eq. |2] that for j = I the final trace over 
the environment system is equal to one and therefore we 
always have fjj = 1. Thus, this decoherence model only 
affects the off-diagonal terms in the <7 z -basis, in other 
words, the oveigenbasis is a pointer basis. 

The remaining task is to evaluate the decoherence fac- 
tor /oi (n, T) since on general grounds one can show that 
fji{n,T) = fij(n,T)*. To evaluate fji(n,T) it is conve- 
nient to notice that the integrals in Eq.|^can be brought 
forward through the independent operator terms in the 
sequence in Eg. 1221 and the evolution can be expressed as 
the succesive application of a superoperator on the initial 
density matrix of the environment p E (Q). Thus, we can 
write: 

f 01 (n,T)=Tr E [O n (p E (0))] (23) 
where the superoperator O is defined as 



O(p) 



f ^K E V E p {V E )\K E V 



-2inu s T/n f j e g v -«r(g+njs)T<r,/r. 



xpe 



2a 



i7r( %-v E )Ta z /n p ierr y 



The dependence of /oi on the self-evolution of the system 
factors out as a phase factor that modulates the overall 
evolution in Eq. This trivial phase factor will be 
omitted from here on because it can be easily restored 
if necessary. After integrating over the random variable 
the last expression becomes 

O(p) = C L-i^+^)T^/n pe -i^-u E )T/n^ + (25) 

where 7 = c — d — sin(2a)/2a and c+ d = 1. It is worth 
stressing that this superoperator is not trace preserving 
or hermitian. It is easy to show that 0(o~ x ) = cr x and 
0(a y ) = j<T y . Thus, both a x and a y are eigenvectors of 
the superoperator O (respectively with eigenvalue 1 and 
7). We will later need to find the other two eigenvectors 
which are linear combinations of the identity (/) and o~ z . 
The decoherence factor /oi is given following a final trace 
over the environment qubit. So, the traceless terms in 
p E (Q), those proportional to a x and a y do not contribute 
to the final result. Thus, to compute foi(n, T) the super- 
operator O is applied n-times to the part of the initial 
state with components along the identity and a z . Writ- 
ing the initial density matrix of the environment qubit as 
p E {0) = (I + p x o- x +PyO- y +p z a z )/2 we obtain 

f 01 (n,T) = \Tr{O n {I)) + l - Pz Tr{O n {a z )). (26) 

The action of O on the identity and a z is 

0(1) = cos(7rOT/n)J-i7sin(7rOT/n)<7 z 
0(<t z ) = ~isin(TrnT/n)I + -fCOs(TTilT/n)a z . 

Note that the above expressions have no dependence on 
the frequencies v$ and v E since they came in as trivial 
phase factors. 

The eigenvalues Ai and A2 (and the corresponding 
eigenvectors) of the superoperator O can be obtained di- 
rectly, giving, 



1 /(I + 7) 2 

Ai = -(l+7)cos(7r»T/n)±W v " cos 2 [ttVLT / n) -7 

(27) 

and, from them one can find the following exact solution, 
cos(7r^T/n)(A5 l - A£) + AiA£ - A 2 A™ 



/oi (n,T) = 



ip z sm(nflT/n) 



(Ai - A 2 ) 
(A? - XS) 



(Ai - A 2 ) 



(28) 



(24) 



Notice that this formula is an explicit expression (ob- 
tained with no approximations) valid for all values of the 
parameters defining our model (n, 7, etc). Also, it is 
worth stressing that the dependence on the initial state 
of the environment (entering the above equation through 
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the initial polarization p z ) is rather trivial. Moreover, 
the first and second lines of the last equation clearly sep- 
arate the real and imaginary parts of the decoherence 
factor /oi. Below, we will analyze the predictions of this 
model for some simple cases. 

Dependence on Kick Angle: Limiting Cases 

We first consider the dependence of the decoherence 
factor /pi on 7. Let us consider three cases. First we 
discuss the limit 7=1 that corresponds to unitary evo- 
lution (that is, no kicks since the kick angle a — 0). Then, 
we consider the case 7 = that corresponds to averaging 
over angles between and 2ir. Finally, we analyze in some 
detail the case where 7 is close to one (small angle kicks), 
which is the condition met in our simulations and experi- 
ments. In all these cases the decoherence factor /01 is di- 
rectly related to observable quantities (<rf ) = 23?[poi/oi] 
and (o£) = 29f[poi/oi]. 

Unitary evolution: 7 = 1 



S 250 - 
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This is the simplest case. Here, the superoperator O is 
such that O(p) = pc?q>(—iirQTo- z /n) for any operator p 
that is a linear combination of the identity and o~ z . (Note 
that we showed earlier that a x and a y are eigenvectors of 
O and thus vanish after the trace). Using this, or simply 
replacing 7 = 1 in the decoherence factor (Eq. 128(1 , 



/oi(n = 0,T) = cos(TrrJT) - i p z sin(vTftT). 



(29) 



FIG. 3: The decay rate as a function of the kick rate F. For 
Q = 300Hz and 7 = 0.98 the kicking is no longer effective at 
inducing decoherence beyond a kick rate of about 50 kicks/ms. 
Only kick rates upto 1000 kicks/ms are shown and after 5000 
kicks/ms the decay is no longer exponential. Inset: The decay 
rate as a function of the kick rate is linear for small values of 
T. The plot is for Q. = 300Hz and 7 = 0.98. 



This has a clear physical interpretation. The decoherence 
factor is independent of the kicking rate (as it should 
be since there are no kicks in this limit). Recall that 
p z is the initial polarization of the environment qubit; 
therefore the system qubit rotates independently of the 
environment qubit. 

Complete randomization: 7 = 

Here the kick angles tj vary over the entire interval be- 
tween and 27r. In this case the above formulae simplify 
substantially to 



/oi(r,T) = cos 1J ( 



■kVL. 



1 p z sm( 



H-W). (30) 



r 



where use of n — FT was made. (Recall that V = n/T 
is the kick rate). In the large T limit we clearly see a 
Zeno-like effect (for an operator not a state) that can be 
obtained from Eq. |30l by noting that 



cos 



K r > 



(1 



1 / n ^\2\rT 

2 F 



, (vrm 2 T, , . 
° XP( ~ 2T } - (31) 



Thus, in this limit for faster kick rates T the system takes 
longer to decohere. 



Average over small angles: 7 = 1 — 0(a 2 ) 

Here we consider the case where the averaging is over 
small angles (the regime we consider in the simulations 
and experiments is a — k/20), where 



7 : 



(32) 



Defining e = 



2 „,2 



we can expand both eigenvalues in 



powers of e to obtain an expression which is valid for 
small n = TT: 



/oi(r,r) = 

(1 - -) rT 
v 2' 



(i 



-)[cos(7rOT) 



i Pz sm(7rOT) + e>(e)]. 



(33) 



In this regime the envelope of the decay of foi is expo- 
nential with a decay rate proportional to the kick rate 
because e -C 1 implies (1 — |)™ w exp(— ne). The analy- 
sis of the exact formula shows that in this case (large n) 
a Zeno-type effect arises (as before). 

The dependence of the decay rate, T2, as a function of 
the kicking rate is shown in Fig. [3] The numerical data 
in the figures are obtained from the exact expression for 
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foi- The initial state for the system qubit is taken to 
be p s = \(I + a x ), in which case /oi is directly propor- 
tional to the transverse polarization of the qubit. For 
small values of the kick rate 1/T 2 is linear in T. However 
for larger values 1/T 2 saturates and decays again due to 
the Zeno-like effect. These results substantiate our ex- 
pectation that the kick-rate can be applied to control 
the attenuation of the recurrences. In the low kick-rate 
limit the role of the kick-rate is analogous to the variable 
environment size in Zurek's model. 

Kraus Forms 

For a one-spin environment a phase dam ping channel 
can be represented by a purification basis |23l |24| that 
evolves the system and environment with the unitary op- 
erator: 



Use = e 



= E s + e- ld <' 2 +E s _e w <' 2 , (34) 



where E^. = \(I ± erf ), or equivalently E + = |0)(0 
and £L = |1)(1|. This operator transforms the states 
of p SE (0) = p s (0) <g> Ef as follows: 



v s = 

vex = 630// z 
ve 2 = —630Hz 
J SEl = 250Hz 
Jse 2 = 50Hz 
JexE 2 = 1747/2 

6 = { — ^jp +;§)} (randomly choosen) 



TABLE I: This table lists the parameters for the model Hamil- 
tonian of Eq. ^3 



the a z a z interaction is sufficient to describe the open- 
system dynamics of phase-damping. This is because de- 
phasing is a special case where the Lie algebra of the 
noise consists of only the two operators a z and / (out of 
a possible four). In contrast, for an arbitrary completely 
positive map the dimension of the environment must be 
at least N 2 for a system with dimension N to induce an 
arbitrary mapping on the system. 



\0)s\0)e cos(0/2)|O)s|O) £ + sm(9/2)\0) s \l) E , 
u 



|l)s|0}i 



cos(0/2)|l) s |O) E - sin(0/2)|l) s |l) B . (35) 



By tracing away the environment states {TLe ■ 
{\0)e, |1)b}) this channel has the Kraus operator sum 
representation 0, 0, |2^| given by 



S{p s ) = MopHln + M lP s M Xl 

where 

M = cos(6>/2)/ 5 , Mi = sin(0/2)of , 
and S is the superoperator map, 

S(P S ) 



Poo PPox 



(36) 



(37) 



(38) 



with (3 = cos 2 (9/2) — sin 2 (9/2). If we parametrize 

cos(0/2) = ^I(i + / 01 ), sin(0/2) = yjl(i-f 01 ), 

(39) 

then we obtain the Kraus operator sum representation 
for the phase damping channel in our model: 

S(P S ) = ~(1 + foi)p S + \(l - kx)a S zP S a s z . (40) 

From the analytical solution to the two-qubit model 
we see that a single qubit environment interacting with 
a single qubit system is sufficient to represent the phase- 
damping channel. Similarly, an iV-dimensional system 
interacting with an environment of dimension TV through 



Two-Qubit Environment: Numerical Simulation 

In the more general case where we wish to implement 
any completely positive map |l8l Il9| on one-qubit sys- 
tem the minimum required environment is two qubits. 
We therefore want to consider a two-qubit environment 
model. Moreover, we want to examine the effect of only a 
finite number of realization of the random kick variables. 
Therefore, a three-qubit model is explored numeric ally 
below. The results of on an NMR QIP simulation [20ll2lT| 
of this model are presented in the next section. 

We now consider the following system-environment 
Hamiltonian, 



where, 



7~Ltot — 7~Ls + He + Use + 'He 1 e 2 i 



7i s — tivs<j z , 

2 



(41) 



H 



E 



7~Lse 
T~Le 1 e 2 



n^VEiO z \ 

i=l 
2 

7T x ^ o 
2 2-^ J SEi<! z O 



Ei 

z ' 



i=l 



E 2 
i ' 



The environment spins, E\ and E2, are also subjected to 
periodic, instantaneous kicks with an evolution operator 
of the form, 



K m = exp 



(42) 
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FIG. 4: Some example decays of the system coherence given 
by {(0|p s |l)} obtained from numerical simulation using Mat- 
lab for 50 realizations of the random kick variables. The kick 
rates for each subplot are labeled above the figures. In this 
range the envelope of the decay is exponential (see inset to 
Fig. |SJ. We note that a higher kick rate leads to a faster 
system decay. 




time(ms) 

FIG. 5: A numerical simulation to demonstrate the suppres- 
sion of revivals at longer times and higher averages. The times 
go out to 500ms and the averages are taken for 200 realiza- 
tions. Note that the revivals that seem prevalent in Figs. [I] 
and |S]are diminished. 



where the 6™ are the random values of the m'th kick. The 
instantaneous nature of the kicks allows the evolution of 
the full system over the time interval T, with n instaneous 
kicks, to be described by the operator, 

U n = fc n exp[—iHtot(T/n)]JCn-iexp[—iHtot(T/n)] 
x ■■• x K ie xp[-fHtot(T/n)], (43) 

where each IC m has a different random kick variable. 

The resultant system density matrix for a single real- 
ization is now obtained by tracing out the environment. 
As before, we are interested in the system coherence as 
expressed through the off-diagonal elements of the sys- 
tem state in the basis of the pointer states, 

(rV(T)|l) = (0\Tr ElE2 [U n p SE ^(0)Ut] \l) a , (44) 

Finally, we must average over different realizations of the 
random variables, which gives the quantity, {(0|p s |l)}, 
where the curly brackets on (0\p s \l) denote the average 
over the finite number of realizations. 

In order to simulate the physical system used in the 
NMR study, we have selected the parameter values pre- 
sented in Table [I] The system and environment were 
initialized in the state <rf E^ 1 and we simulated the 
evolution of the system on Matlab. We ran 10 different 
kick rates that ranged from 3 kicks/ms to 30 kicks/ms in 
steps of 3. The kicks were sampled from a uniform dis- 
tribution of angles that ranged between — 5j to ^tj. The 



series was run for 150 ms. We averaged over 50 realiza- 
tions and obtain the plots shown in Fig. 0] As shown in 
Fig- El the late-time oscillations reflect the finite number 
of realizations of the random variables. The envelope 
of the decays in Fig. 21 were fit to an exponential and 
the decay constants exhibited a linear dependence on the 
kick rate for small kick rates, as expected from the an- 
alytic solution (see Fig. EJ. At about 900 kicks/ms the 
decay rates start decreasing with increasing kick rate and 
the system starts to become decoupled from the environ- 
ment, an effect noted earlier in Eq. 1301 This is the well- 
known decoupling phenomena in NMR [23. The onset 
of decoupling occurs when the rms angle of the stochas- 
tic kicks approaches a rotation of n (criticial damping). 
The rms angle is given by the typical kick size ~ 7r/10 
times the square root of the number of kicks over a cycle 
time ~ 1/2 J of the system-environment interaction. For 
the strongest system-environment coupling, J ~ 250 Hz, 
the onset of decoupling is expected at a kick rate of 800 
kHz, in good agreement with the numerical results (see 
Fig.©. 



IV. THE NMR IMPLEMENTATION 

In this section, we describe the experimental imple- 
mentation of our model. We chose propyne as the physi- 
cal system (see Fig.Qfor the internal Hamiltonian param- 
eters). The hydrogen indicated with a circled 1 represents 



9 



350 



300 



250 



100 



30 








* 25 
20 










E 
t 

10 






» 5 
# 






» 5 10 15 20 25 30 
I Number of Kicks/ms (r) 

Til 



Number of Kicks/ms (r) 



FIG. 6: Numerical simulation of the decoherence rate and 
the decoupling limit. Beyond a kick frequency of about 900 
kicks/ms the decoherence rate from the kicking starts to de- 
crease. This transition to a decoupling effect is described in 
the text. After about 5000 kicks/ms the decays are no longer 
exponential. Inset: Demsontration of the proportionality be- 
tween decoherence rate and kick rate for low kick rates. This 
linear relationship can be understood from the analytic results 
obtained for the one-qubit environment. 



the system qubit and the two carbons labeled with a cir- 
cled 2 and a circled 3 represent the environment qubits E\ 
and E2, respectively. These spin-1/2 nuclei have a large 
resonance frequency offset, so the hydrogen and carbon 
can be addressed and detected separately. The relatively 
large couplings present amongst these nuclei implies the 
interactions take place over short times, and the long re- 
laxation times allow one to observe the hydrogen signals 
over a relatively long time span without significant nat- 
ural decay. The experiments were carried out on a liquid 
solution of propyne using a Bruker Avance spectrometer. 
Neglecting the methyl group (because it couples in very 
weakly) , the internal Hamiltonian for propyne is given to 
a good approximation by 

Hint = t^[v\o\ + viol + v z a\ 

+|(Ji2<^ + J 23 cr 2 • o z + J 13 ala 3 z )}, (45) 

where the v y s are Larmor frequencies and the J's the 
spin-spin coupling constants in Hertz (the various values 
are given in Fig. EJ). Eq. 021 should be compared with 
Eq. 1411 The non-secular coupling between the carbon 
nuclei can be observed in the carbon spectra but has a 



H 




FIG. 7: The propyne molecule. The encircled labels on the 
13 Cs and the rightmost hydrogen index the spins used in the 
experiment. The methyl group consists of the three hydrogens 
and an unlabeled carbon. In the experiments the field of 
the spectrometer was ~ 9.2 T and the hydrogen resonances 
were ~ 400 MHz while the carbon resonances were around 
~ 100 MHz. The chemical shift difference between the two 
labeled carbons is 1.260 kHz. Using the indexing scheme in 
the figure the J-coupling constants are as follows: J12 = 246.5 
Hz, J 23 = 173.8 Hz and J 13 = 51.8 Hz. The longitudinal 
relaxation times are T{ = 8.7s, If = 23s and Tj 3 = 43s, while 
the transverse relaxation times are T\ = 1.1s, T-f = 1.9s and 
Ti = 1.7s. 



negligible effect on the relevant experimental results. 

A convenient choice for the initial state of system and 
environment is one where hydrogen is in a superposition 
state and both carbons are in an eigenstate. By placing 
the methyl hydrogens in an eigenstate as well, they can 
be eliminated from playing a role in the hydrogen spin 
dynamics. This was accomplished by using a highly selec- 
tive rf pulse that irradiated a spectral line corresponding 
to the state 

a^E^E^E™, (46) 

where E + = \{I + <J Z ), H represents hydrogen, CI car- 
bon 1, C2 carbon 2 and M the methyl hydrogens. For 
this implementation we used a 5.5s EBURP1 [13, Ei| 
pulse. The spectral resolution of this pulse was .5 Hz 
and its design is such that it only generates a uniform 
excitation profile in the specified bandwidth. Ultimately, 
only ~ l/10th of the maximum intensity was excited. 
Nonetheless, this yielded sufficient signal-to-noise ratio 
to carry out the experiments. 

The observed hydrogen signal corresponds to (a^ (t) + 
i(jy(t)), and is equivalent to tracing away the carbons. 
The peaks of hydrogen spectrum had linewidths of ~ .4 
Hz. Consequently, the hydrogen signals decayed very 
slowly and we were able to pick a 150 ms portion of the 
absolute magnitude that remained flat within one per- 
cent. 
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FIG. 8: Example decays from the experiment. The hydrogen 
signal was directly detected and the real part of the complex 
signal is plotted. The fluctuations at the tail end of the higher 
kick rates are due to low statistics. This was confirmed by 
comparing with simulations at higher averages. (See Fig. |5J 
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FIG. 9: The linear dependence of the experimental decay 
constants on kick rate. The data point symbols (x) are larger 
than the error bars which range from ±.0075 to ±.0573. Com- 
pared to the slope in the simulations of Fig.|S|the experimental 
slope reflects faster decay. This disparity is due to the slight 
differences between the experiment and simulations. 



The carbon spin dynamics consisted of a series of 
delays interleaved with pulses. During the delays the 
spins evolved under the internal Hamiltonian. The pulse 
flip angles were randomly sampled from a uniform dis- 
tribution that ranged between — 7r/20 to 7r/20 about 
the y— axis. A cycletime of 1 ms was defined within 
which the kick frequency ranged from 3 kicks/cycle to 
30 kicks/cycle in steps of 3 for a total of 10 different kick 
frequencies. The range of the kick frequency was limited 
by the shortest pulse the spectrometer was capable of 
generating, which is 100 ns. The time alloted for a se- 
quence of one delay period followed by a pulse was given 
by the cycletime/ (number of kicks/cycle). Within this 
sequence the delay time is given by the total sequence 
time minus the pulse-on time. The maximum pulse-on 
time was 10 jis which corresponded to the maximum flip 
angle of 7r/20. The nutation frequency for this RF field 
was 2500Hz. (Compare this to the chemical shifts of the 
carbons which were separated by 1260Hz). For a given 
kick frequency, the length of the series of successive se- 
quences of delay plus pulse, generated as described above, 
fit the total acquistion time of 150 ms. 

The experiments were run for 10 different kick frequen- 
cies with an average over 50 realizations. A waiting time 
of 300s was used between successive realizations. Fig. |H] 
shows the result of the experiments. The absolute mag- 
nitude of these plots were fit to an exponential. The x 2 
per degree of freedom for the 10 fits ranged from 1.1 to 
8.5. The x 2 s fit the average decays well but don't ac- 



count for the details in the fine structure evident from 
the oscillations of the magnitudes. As the kick frequency 
increases the data demonstrates that the system is deco- 
hered faster. A plot of the decay constants as a function 
of kick frequency, Fig. [§1 shows this trend clearly. The 
experiment results seem to exhibit revivals in the higher 
kick rates of Fig. But this is due to low statistics (see 
Fig. E| . 



V. DISCUSSION 

We have described a method for modeling decoherence 
that requires only limited quantum resources, and im- 
plemented the model on an NMR QIP. The key feature 
of the model which enables simulation of the dephasing 
effects and the attenutation of recurrences normally pro- 
duced by a much larger quantum environment is the ap- 
plication of classical kicks to randomize the information 
in the environment states. Although the quantum sys- 
tem and environment dimensions are small and remain 
fixed, the system state exhibits an irreversible loss of co- 
herence due to an averaging over the random realizations 
of kicks to the environment states. In particular, in the 
case of a <r z a z system-environment interaction we have 
shown that the kick frequency can be varied to control 
the decay rate of the phase-damping. Although in this 
paper we have focussed on the simulation of continuous 
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phase-damping, the model can be immediately general- 
ized to other system-environment couplings and the re- 
sultant decoherence channels. A major advantage of this 
model is that it provides a procedure through which the 
mechanisms of decoherence can be explored using tech- 
niques currently available in NMR QIP. 

As resources permit, the model we have described may 
be generalized to simulate and study a wider variety of 
decoherence channels and system-environment couplings. 
In particular, the "nearest" quantum environment need 
not be the only quantum environment. For example, in 
order to implement a time-varying decoherence process 
with a fixed set of system-environment couplings it may 
be advantageous to introduce an environment "hierar- 
chy" (see Fig. 1 for a schematic). The idea here is to 
couple the first quantum environment to a second, larger 
environment (through another set of fixed bilinear cou- 
plings), and so on. The dimension of the next Hilbert 
space in the environment hierarchy may be limited to 
N 2 , where N is the dimension of the Hilbert space of the 
immediately smaller system. In this framework only the 
nearest environment remains directly coupled to the sys- 
tem of interest. The approximation of using stochastic 
classical fields to reduce unwanted back-action may then 
be applied to the final quantum environment, which is 
much more remote from the system of interest. 

In conclusion, we have developed a model that is prac- 
tical for simulating quantum decoherence effects associ- 
ated with a time-independent superoperator on a QIP de- 
vice. By varying the phase kicking rate in the stochastic 
Hamiltonian we can control the system's phase-damping 
rate. In this presentation we have shown the effectiveness 
of the methodolgy, in the case of one and two spin en- 
vironments, using analytical solutions, numerical simula- 
tions, and a physical implementation on an NMR QIP de- 
vice. These methods have illustrated the use of stochas- 
tic kick rates on the quantum environment for controlling 
system decoherence rates and recurrence times. 
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